511 keV Photons from Superconducting Cosmic Strings 
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We show that a tangle of light superconducting strings in the Milky Way could be the source 
of the observed 511 keV emission from electron-positron annihilation in the Galactic bulge. The 
scenario predicts a flux that is in agreement with observations if the strings are at the ~ 1 TeV 
scale making the particle physics within reach of planned accelerator experiments. The emission 
is directly proportional to the galactic magnetic field and future observations should be able to 
differentiate the superconducting string scenario from other proposals. 
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The detection of a bright 511 keV line by the SPI spec- 
trometer on the International Gamma-Ray Astrophysics 
Laboratory (INTEGRAL) satellite, has established the 
presence of a diffuse source of positrons in the Galactic 
Center (GC) Q. The observed photon flux of 
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with a linewidth of about 3 keV is in good agreement 
with previous measurements Q. The observations sug- 
gest a spherically symmetric distribution. Assuming a 
Gaussian spatial distribution for the flux, a full width at 
half maximum of 9° is indicated. 

The origin of these Galactic positrons remains a mys- 
tery. Several scenarios involving astrophysical sources 
have been proposed, including neutron stars or black 
holes, massive stars, supernovae, hypernovae, gamma- 
ray bursts or cosmic rays |3|. However, the fraction of 
positrons produced in such processes is uncertain, and it 
is unclear that the positrons could fill the whole bulge. 

Alternatively, mechanisms associated with the Dark 
Matter (DM) at the GC have been put forward. If DM 
is constituted by a light (1-100 MeV) scalar, its anni- 
hilations could account for the observed signal Q. An- 
other possibility is that DM could be in the form of non- 
hadronic color superconducting droplets. Positrons on 
the surface of the droplet could annihilate with ambient 
electrons producing a gamma-ray line [5j . 

In this Letter, we show that a network of light su- 
perconducting strings [6j occurring in particle physics 
just beyond the standard model could also be a source 
of positrons. More detailed observations of the positron 
distribution could be used to distinguish this source from 
the other possibilities. Moreover, this resolution of the 
positron observation works with strings at an energy scale 
of ~ 1 TeV, since the flux from heavier strings is lower. 
Therefore, experiments at the Large Hadron Collider in 
the near future will also probe particle physics relevant 
for this explanation of the positrons. 

We assume that a tangle of light superconducting 
strings exists in the Milky Way. The strings are frozen in 
the Milky Way plasma as long as the radius of curvature 
is larger than a certain critical length scale. If the cur- 
vature radius is smaller, the string tension wins over the 



plasma forces and the string moves with respect to the 
magnetized plasma. During the string motion, the loop 
will cut across the Milky Way magnetic field, generating 
current as given by Faraday's law of induction. The cur- 
rent is composed of zero modes of charged particles, in- 
cluding positrons, propagating along the string. Once the 
energy of a positron on the string exceeds the rest energy 
of 511 keV, it becomes possible for it to leave the string 
primarily due to scattering by counter-propagating parti- 
cles. Even though the scattering rates are model depen- 
dent, general arguments discussed below indicate that 
the scattering of positrons with counter-propagating light 
quarks will be very efficient at ejecting the positrons at 
the threshold of 511 keV. The ejected positrons will anni- 
hilate with the ambient electrons, thus emitting 511 keV 
gamma rays. 

The string dynamics is determined by comparing the 
force due to string tension to the plasma drag force. The 
force per unit length due to string tension fi is 
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where R is the string radius of curvature. The drag force 
is more complicated, needing an analysis of the plasma 
flow around a current carrying string 0, . The string 
effectively behaves like a cylindrical body around which 
the plasma flows. The radius of the cylinder, in nat- 
ural units, is r s = J / (v rc \y/p) , J being the current on 
the string, v rc i the velocity of the string relative to the 
plasma, and p the plasma density. The drag force on 
such a cylinder is 



F d ~ pv^ cl r s ~ yfpv Te iJ. 
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The maximum value of the damping force occurs when 
v rel = 1. Therefore, F d>max ~ JpJ. If F s > Fd, maX) the 
string moves under its own tension and the damping force 
is insufficient to prevent the string from accelerating to 
relativistic velocities. String loops will then emit electro- 
magnetic radiation and eventually dissipate. If, however, 
F s < -Fd.max, the string will accelerate until the damping 
force grows so as to cancel the tension. Then the string 
moves at a terminal velocity relative to the plasma found 
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by equating F s and F^: 
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Setting Wtcrm = 1, we find the critical radius of curvature 
R c when the strings move relativistically: 
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To summarize, strings with tension p, current J, and 
curvature radius R > R c , move at v tC i m with respect to 
a plasma of density p. If R < R c , the string moves at 
relativistic speeds. 

In a turbulent plasma, such as in our Milky Way, there 
is another length scale of interest, called i?„ > R c ), 
even when the string motion is overdamped. For R > i?*, 
the terminal speed of the strings is small compared to the 
turbulence speed of the plasma and the strings are car- 
ried along with the plasma. As the strings follow the 
plasma flow, they get more entangled due to turbulent 
eddies, and the strings get more curved until the curva- 
ture radius drops below Then the string velocity is 
large compared to the plasma velocity, and, hence, the 
strings break away from the turbulent flow. Therefore, 
i?* is the smallest scale at which the string network fol- 
lows the plasma flow. For i?„ > R > R c , the string 
motion is over-damped but independent of the turbulent 
flow. Hence, string curvature on these scales is not gen- 
erated by the turbulence, and we can estimate the length 
density of strings in the plasma as pi ~ To find 

i?*, denote the turbulent velocity at the scale i?„ by v* 
and set wtcrm ~ w*, to get \u\ 
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where, for convenience, the dimensionless parameter k 
has been introduced via J = Ke^/p. 

To determine v*, we note that, on scales less than 
I r~j lOOpc = 3 x 10 20 cm the interstellar plasma has 
velocity spectrum given by magnetohydrodynamic turbu- 
lence as vr ~ ^/(i?//) 1 / 4 , where v\ ~ 10 6 cm/s. Inserting 
this expression in Eq. (Jfjj and solving for i?* gives 
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Note that in these estimates we have assumed R* < I; 
otherwise, the turbulent scaling law will be different. 

The current on the string is generated by the motion of 
the string across the galactic magnetic field as described 
by Faraday's law of induction and is limited by the mi- 
croscopic interactions of the charge carriers on the string. 
A given charge carrier can, in principle, leave the string 
once it has enough energy. The escape is triggered by 
several factors, of which the most important is scattering 





off counter-propagating zero modes (bosons or fermions) . 
The current in any particular species saturates once the 
scattering rate equals the growth rate due to Faraday in- 
duction. The rate of current growth for particle species 
(denoted X) with charge e is given by Jx ~ e 2 vB. The 
time scale for current growth, assuming v ~ c, is 



growth 
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The current decays because the charge carriers on the 
string can interact and leave the string. The positron 
current is limited mainly by scattering with counter- 
propagating quarks {u quarks for electroweak strings). 
Since the lightest quarks have ~ 1 MeV mass, which is 
similar to the positron mass m e , the decay time when 
the positron current is order em e is r e = 1/na, where 
n ~ (J e /e)m 2 is the number density of positrons and 
a ~ a 2 /m 2 is the electromagnetic scattering cross sec- 
tion, a being the fine structure constant. Therefore, 
r e ~ ea~ 2 J^ 1 ■ Equating r e and the growth rate given by 
growth. ^ get J e ~VB~ 10- 4 eV, with B ~ KT 6 G 
and k ~ 1. Since this is much less than m e = 511 keV 
and positrons cannot escape the string until they reach 
511 keV energy, the threshold for their escape is also 511 
keV. Here we should note that the quark that scatters 
off the positron cannot completely escape from the string 
unless it gains ~ 100 MeV energy and hadronizes. How- 
ever, at energies less than 100 MeV, the quark can still 
get kicked out from its zero mode state on the string to 
a quark that is not in a zero mode. Such a quark will be 
confined to a (Aqcd) -1 ~ (100 MeV) -1 shell around the 
string. Pion emission from the string cannot deplete the 
baryonic current on the string. Only at ~ 1 GeV energies 
can antiprotons be emitted, leading to anoth er p ossible 
signature of galactic superconducting strings [l^] • 

As a piece of string of length R„ cuts through a mag- 
netic field B, it will produce electrons or positrons, with 
equal likelihood, at the rate dN/dt ~ ev+BR*. In a vol- 
ume V = 4tt j L 3 /3, there are - L 3 /Rl such pieces of 
string, and, hence, the rate of particle production in the 
entire volume is 
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The current in the positrons will grow at first but then 
saturate at 511 keV. After that, further motion of the 
string across the galactic magnetic field will generate 
positrons that leave the string. So Ny is also the number 
of positrons being produced in the volume V which we 
denote by N + . Inserting Eq. in |JSJ we get 



dN^ 
dt 
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We are interested in a region of radius 1 kpc around the 
galactic center and so L 3 ~ Vi(lkpc) 3 , where V\ is a di- 
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mensionless parameter. The plasma density in the galac- 
tic center is higher than the average Milky Way den- 
sity. For an estimate, we adopt an isothermal sphere 
scaling normalized to p§ ~ 10 _25 gm/cm 3 at 8 kpc from 
the center, thus giving p ~ 6 • 10~ 24 p gc gm/cm 3 , where 
we have introduced the parameter p gc . The magnetic 
field in the galactic center is not known very precisely 
but is estimated to be ~ 10 -3 G Therefore, we 

will set B = i?3l0~ 3 G. (Note the conversion of mag- 
netic field strength in Gauss to particle physics units: 
1G = 1.95xl0" 20 GeV 2 .) The parameters v t and I will be 
different in the galactic center but these are not known. 
We will write vi = 10 6 vi$ cm/s and I = 100 Zioo P c - The 
string tension will be taken to be fi = TeV) 2 , where 
Hi is a parameter, and e 2 ~ 0.1, and then Eq. (fTO)) gives 
us 

dN+ 



dt 
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Although the astrophysical parameters describing the 
galactic center are not known very accurately, assum- 
ing equipartition of plasma kinetic energy (~ pvf) and 
magnetic energy (~ B 2 /8ir), with I ~ p -1 ^ 3 , we find that 
vifi ~ 100 and Zioo ~ 0.1, which boosts the estimate in 
Eq. HIl) by 10 5 , yielding 

^ < 10 47 s" 1 . (12) 
dt 

From Eq. J7J), the string velocity with respect to the 
plasma is v* ~ 10~ 5 c. Hence, the positrons will leave 
the string with a small Lorentz factor, 7 ~ 1, and they 
will gradually slow down by Coulomb collisions in the 
interstellar medium (ISM). The energy loss rate is ap- 
proximately [3 dE/dt - 2 x 10" 9 (iV ff /10 5 m" 3 )(log7-|- 
6.6)eV/s, where Njj is the number density of target 
atoms. This rate yields a stopping distance of 10 24 cm. 
Considering a simple random walk, with the Larmor ra- 
dius on the order of 10 1 cm, the positrons travel a distance 
~ 10 12 cm, so they are easily confined to the galactic cen- 
ter. Note that this distance is less than the typical string 
separation estimated from Eq. J7|), R* ~ 5 • 10 15 cm, so 
the positrons will have linclike features on angular scales 
~ 10~ 7 , which is too small to be resolved with the 2° 
angular precision of INTEGRAL. Once produced, the 
positrons will undergo different processes in the ISM. Pair 
annihilation with an ambient electron and positronium 
(Ps) formation and decay via para-Ps will occur both in- 
flight and after thermalization. The resulting spectrum 
depends on the specific details of the ISM, but detailed 
analysis shows that a narrow 511 keV line generically re- 
sults, AE ~3 keV, in agreement with observations [l5| . 
Hence, the spectral shape does not depend on the details 
of the strings, and in this respect the model cannot be 
distinguished from alternative mechanisms such as light 
DM or astrophysical sources. 

Three-quarters of ortho-Ps annihilate in a 3 photon 
continuum final state. As a result, each positron will 



contribute (2 — 3/p s /2) 1 511 keV photons, where the 
Ps fraction has been measured to be fp s = 0.93 ± 0.04 
for the galactic center 16]. Multiplying the observed 
gamma-ray flux given in Eq. by the area of the sphere 
at our location (8 kpc) around the galactic center, we get 
the actual positron production rate in the galactic center: 
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dt 
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1.2 x 10 43 s 
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Comparing Eqs. (|12H and (|13fl . we conclude that light 
superconducting strings are possible sources of positrons 
that lead to the flux of 511 keV gamma rays observed by 
the INTEGRAL collaboration. 

We see from Eq. (|11|) that a unique prediction of our 
scenario is that the gamma-ray flux is proportional to the 
magnetic field strength in the Milky Way, with a milder 
dependence on the plasma density. In the disk, the mag- 
netic field intensity decreases by B3 ~ 10~ 3 . At the 
same time, looking towards the disk, the volume of the 
Milky Way that contributes positrons is larger. Taking 
these factors into account, for a disk with thickness 1 kpc 
and a radial extent of 30 kpc, we estimate a photon flux 
~ 10~ 6 cm _2 s _1 in a 16° field of view as in SPI. In a di- 
rection perpendicular to the disk, the volume contribut- 
ing to the flux will be smaller, adding to the suppression 
of the magnetic field to yield a flux of ~ 10 -7 cm _2 s _1 . 
Thus far, there have been no reliable detections outside 
of the central region of our Galaxy, with SPI placin g an 
upper bound of 1.2 x 10~ 4 cm~ 2 s _1 on the flux |17| . 
The target sensitivity of the SPI instrument, once suffi- 
cient exposure becomes available, is 2 x 10~ 5 cm~ 2 s at 
511 keV 0], somewhat above what is needed to map 
the emission from the disk in our scenario. 

That the flux should follow the magnetic field is in 
marked contrast with the MeV DM hypothesis. There 
the flux follows pp M , and a signal from nearby DM dom- 
inated regions, e.g., the Sagittarius dSph galaxy, is ex- 
pected [l!| . No magnetic fields have been measured 
in Sagittarius, although low surface brightness galaxies, 
which are somewhat similar, show \xG scale magnetic 
fields [iof. Therefore, if superconducting strings source 
the observed 511 keV, the estimate would be weaker by 
10 -3 than that of Eq. I|ll|) due to the weaker magnetic 
field, and another factor ~ 10 _1 due to the source being 
3 times further away. The volume of Sagittarius dSph 
is comparable to the Milky Way bulge, leading to, at 
most, a flux of ~ 10~ 7 cm~ 2 s _1 in the direction of Sagit- 
tarius, some three orders of magnitude fainter than the 
MeV DM model prediction. After two Galactic Center 
Deep Exposures, INTEGRAL has not detected emission 
from Sagittarius, although the effective observation time 
is not yet sufficient to reach the sensitivity of the MeV 
DM predicted fluxes [2l| ■ For the color superconducting 
DM scenario, the flux follows p v isibiePDM which is also 
different from our proposed scenario. 
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Present-day galactic magnetic fields could have been 
produced by amplification of a tiny primordial field by 
a galactic dynamo. In that case, currents could build 
up in the network of superconducting strings at earlier 
epochs. Since there is no turbulence at last scattering, 
a tangle of strings will not be formed; there is only ~ I 
string per horizon which would simply be dragged by 
the plasma. During recombination, the plasma density 
drops, and the damping force on strings is vastly re- 
duced. String tension causes strings to move at relativis- 
tic speeds cutting across magnetic fields and generating 
a flux of positrons, which will once again annihilate with 
ambient electrons to produce 511 keV gamma rays. By 
the present epoch, the gamma rays would have redshifted 
to 511 eV lines. One string per horizon gives a positron 
flux of ~ 10~ 17 cm~ 2 s _1 sr , much smaller than the up- 
per limit on the diffuse extragalactic component by the 
Chandra experiment of 10 _3 s _1 cm~ 2 sr _1 at 511 eV, 
assuming a redshifted width of 3 eV. 

Recently, the High-Energy Antimatter Telescope 
(HEAT) balloon experiment has confirmed an excess of 
high-energy positrons in cosmic rays at energies around 
10 GeV |23j. Direct production of 10 GeV positrons from 
superconducting strings does not seem likely, because it 
is hard for the positron current to build up beyond 1 MeV 
in the presence of counterpropagating zero modes. How- 
ever, superconducting strings at the TeV scale could still 
produce heavy charged fermions, if such zero modes exist, 
which could then decay to give positrons in the 10 GeV 
energy range. The production rate of heavy fermions can 
be computed from Eq. lfTT|l to be ~ 10~ 23 cm -3 s _1 . Al- 
though the actual positron yield depends on the branch- 
ing ratio in the particle physics model and the details of 
the diffusion process, comparing this value to the typical 
values of ~ 10 -27 cm _3 s _1 for positrons from annihila- 
tions of dark matter particles at the electroweak scale |24j 
shows that this possibility deserves further investigation. 

We conclude that light superconducting strings could 
produce enough positrons in the galactic center to ex- 
plain the flux of 511 keV gamma rays observed by the 
INTEGRAL collaboration. The scenario can be differ- 
entiated from other proposals by higher resolution ob- 
servations and by observations in directions away from 
the galactic center. Light superconducting strings might 
also produce stellar and cosmological signatures. Since 
the strings may be at the 1 TeV energy scale, the involved 
particle physics is within the energy range of planned ac- 
celerator experiments. 
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